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Abstract 

We study the cosmological evolution of a type-0 string theory by employing non- 
criticality, which may be induced by fluctuations of the D3 brane worlds. We check the 
consistency of the approach to 0(a') in the corresponding cr-model. The ten-dimensional 
theory is reduced to an effective four-dimensional model, with only time dependent fields. 
We show that the four- dimensional universe has an inflationary phase and graceful exit 
from it, while the other extra dimensions are stabilized to a constant value, with the fifth 
dimension much larger than the others. We pay particular attention to demonstrating the 
role of tachyonic matter in inducing these features. The Universe asymptotes, for large 
times, to a non-accelerating linearly-expanding Universe with a time- dependent dilaton 
and a relaxing to zero vacuum energy a la quintessence. 



1. Introduction 

Interest in type-0 string theories has recently arisen due to their connection to four- 
dimensional SU(N) gauge theory via appropriate dualities [|TJ. The type-0 string model 
contains electrically-charged three branes, N of which may be stuck on top of each other to 
construct the dual of a gauge theory. An interesting feature of these non-supersymmetric 
theories is the appearance of tachyonic matter excitations. These, however, do not lead 
to an instability, since they couple to the appropriate five-form field strength that drives 
the tachyon to a non-trivial vacuum expectation value capable of reversing the original 
(negative) sign of the mass-squared. 

Irrespective of the connection with dual gauge theories, the presence of tachyonic matter 
in type-0 string theories has important consequences in implying a cosmological evolution 
of the brane of inflationary type || . The role of tachyonic matter in inducing non-trivial 
cosmological evolution in two-dimensional string models has been studied previously ||, 
with interesting conclusions on the possibility of inflationary phase as well as 'graceful' exit 
from it 0]. It should be stressed, however, that the graceful exit issue has been resolved only 
upon invoking non-criticality of the underlying string theory, which is arguably consistent 
with the non-equilibrium nature of the inflationary phase. An important issue in such non- 
critical string analyses is the identification of target time with the Liouville mode || , the 
consistency of which with the presence of an inflationary phase and the eventual graceful 
exit is a highly non-trivial issue. 

Such a consistency is equivalent to demonstrating the existence of a solution of a system 
of equations which represent a generalization of the conformal invariance conditions in the 
case of a non-critical Liouville string, as a result of the restoration of conformal invariance 
after Liouville dressing ||. An interesting physical issue is what causes the non-criticality of 
the type-0 string. At present, the best scenario we can offer are the quantum fluctuations of 
the three brane solitonic structures 0, |J. In General Relativity there are no Rigid Objects, 



due to general covariance, and so recoil fluctuations of the brane worlds are expected in 
general. Moreover reconciliation with quantum mechanics requires the incorporation of 
position and momentum uncertainties of the brane worlds. 

Under the identification of the Liouville mode with the target time, the induced non- 
criticality leads to interesting predictions on the possibility of an inflationary phase of the 
fluctuating brane ||. However, the interesting issue is the exit from such a phase. It 
appears that, apart from the graviton and dilaton fields on the brane, the issue of graceful 
exit from the inflationary phase requires the presence of additional matter excitations. It 
is the purpose of this work to analyse an effective theory on a three-brane in the context 
of a non-critical type-0 string theory, with the aim of demonstrating such a graceful exit 
from the inflationary phase. As mentioned above, the non-criticality will be associated 
with quantum fluctuations of the embedded three-branes. 

We consider a non-critical ten- dimensional string theory of type-0, and we dimensionally 
reduce the ten-dimensional 0{a!) string-effective action to a four- dimensional space-time 
on a 3-brane world, assuming that all fields depend only upon time. We choose two different 
fields to parametrize the internal space. The first field sets the scale of the fifth dimension, 
while the other parametrize a flat five dimensional space. In the same way, we dimension- 
ally reduce the modified (because of the off-criticality) /3-functions of the ten-dimensional 
theory, to the effective four- dimensional /3-functions, assuming a Robertson- Walker form 
of the four- dimensional metric. Liouville dressing of the non-critical theory should restore 
conformal invariance of the cr-model, which results in a set of generalized conformal in- 
variance conditions. Upon the identification of the Liouville mode with the target time, 
which eliminates scenaria with two target time-like coordinates, such conditions become 
equivalent to 'equations of motion' for the appropriate fields. 

These equations, supplemented by the Curci-Paffuti equation, stemming from world- 
sheet renormalizability of the cr-model, are studied in detail. At present, they appear to 



be too complicated to be solved analytically. We thus follow a systematic method, that of 
quasi-linear systems, to solve these equations numerically. 

The main results of such an analysis may then be summarized as follows: our model 
clearly exhibits an inflationary phase and, most importantly, a graceful exit from it. The 
supercriticality, as well as the tachyonic matter, play a crucial role to this effect. Moreover, 
there is stabilization of the internal space. The fifth dimension contracts during inflation 
and then freezes to a constant value, while the other five internal dimensions also freeze 
to a constant value, which in the solution we present in section four, is different from the 
first, indicating the possibility of obtaining different scales for the extra dimensions. 

The structure of the article is as follows: in section two we give a brief review of the 
Liouville formalism under the identification of the Liouville mode with the target time. In 
section three we discuss the type-0 string theory and we derive the equations of motion of 
the effective non-critical four-dimensional theory. In section four, we describe the method 
for solving numerically the system of the differential equations describing the equations of 
motion and we support our findings by analytic arguments on the asymptotic behaviour of 
the solutions. We present the numerical results and we discuss their physical significance. 
Finally, in section five, we present our conclusions and outlook. Some important formal 
properties of the dilaton field in the context of Liouville strings are discussed in an appendix. 

2. Liouville String Formalism 

We commence our analysis with a brief review of the Liouville dressing procedure for 
non-critical strings, with the Liouville mode viewed as a local world-sheet renormalization 
group scale j5[. Consider a conformal cr-model, described by an action S* on the world- 
sheet E, which is deformed by (non conformal) deformations J s g l Vid 2 a, with V{ appropriate 



vertex operators. 

S g = S*+ \ g % V % d 2 o (2.i; 



s 



The non-conformal nature of the couplings g l implies that their (flat) world sheet renor- 
malization group (3- functions, (3 l , are non- vanishing. The generic structure of such (3- 
functions, close to a fixed point, {g l = 0} reads: 

F = (hi-2)g i + <* jk g>g k + otf). (2.2) 

In the context of Liouville strings, world-sheet gravitational dressing is required. The 
"gravitationally" -dressed couplings, X l (g,t), which from our point of view correspond to 
renormalized couplings in a curved space, read to 0(g 2 ) R [U 



\'(g, t) = g*e^ + Q^fiktfg*^ + 0(g 3 ), Q 2 = - (c - c*) (2.3) 

where t is the (zero mode) of the Liouville mode, c* is a critical (fixed point) value of the 
central charge, about which we perturb the theory, Q 2 is the central charge deficit, and a, 
are the gravitational anomalous dimensions: 

«.(«. + Q) = ^ - 2 for c>c* (2.4) 

Below we shall concentrate exclusively to the supercritical string case, Q 2 > 0, which from 
the point of view of identifying the Liouville mode with target time, corresponds to a 
Minkowskian signature spacetime manifold. The fixed point value c* is normalized such 
that for superstrings we are dealing with here, one has Q 2 = \{d — 9), in the special case 
of superstrings with non critical space time dimensionality (d). 



Due to the renormalization ( |2.3| ), the critical-string conformal invariance conditions, 
amounting to the vanishing of flat-space /3-functions, are now substituted by: 

A* + QX = -/3*(A) for c > c*. (2.5) 

where the notation f3 l {\) denotes flat- world-sheet (3- functions but with the formal substi- 
tution g % — * X l (g,t). Note the minus sign in front of the flat world sheet /^-functions (3 l in 
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(2.5), which is characteristic of the supercriticality of the string || [10|]- As we see later 



the sign will be crucial for the existence of acceptable inflationary solutions demonstrating 
graceful exit from the exponential expansion phase. Upon the identification of the Liouville 
mode with the target time the dot denotes temporal derivative. 

A highly non-trivial feature of the f3 l functions is the fact that they are expressed as 
gradient flows in theory space [O, [EJ, i.e. there exists a 'flow' function F[g] such that 



P = gH 5 Z\A (2.6) 

6g3 



where Q %3 is the inverse of the Zamolodchikov metric in theory space [11], which is given 



by appropriate two-point correlation functions between vertex operators V 1 . In the case of 



stringy cr-models the flow function T may be identified [12j with the running central charge 
deficit Q 2 = c[g] — c*, where c[g] is the central charge of the deformed theory, and c* one of 
its critical values (conformal point) about which the theory is perturbed by means of the 
operators V\ 

It is then interesting to notice that, upon the identification of the Liouville mode with 
the target time ||, the form of equations (|2.5|) , upon the condition (|2.6|) , is reminiscent 
of inflaton equations in standard Cosmology, with the role of the inflaton potential being 
played by the flow function Q 2 [g] in theory space, and that of the Hubble parameter H 
by its square root Q\g\. There is an important difference, though, as compared with the 
inflaton case, in the sense that here all the field modes seem to exhibit this behaviour. 

An important comment we would like to make concerns the possibility of deriving 
the set of equations ( |2.5|) from a target space action. This issue has been discussed in 
the affirmative in fll3| , where it was shown that the set of equations (|2.5|) satisfies the 
Helmholtz conditions for the existence of an action in the 'space of couplings' {g 1 } of the 
non-critical string. The property ( f2.6| ) is crucial to this effect. Upon the identification 
of target time with the Liouville mode || this action becomes identical with the target 
space action describing the off-shell dynamics of the Liouville string. We should stress the 
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fact that the action is off shell, in the sense that the on-shell conditions correspond to the 
vanishing of the /3-functions j3 l . In our case j3 l 7^ 0, and the identification of the Liouville 
mode with the target time implies that the space time graviton /3-function on the right 
hand side of ( |2.5|) , as well as other target-space tensorial structures, appearing inside the 
f3 l functions for the various modes, contain temporal (Liouville) components as well. 

In this respect, our non-equilibrium Liouville string approach to the temporal evolution 
should be contrasted with the interpretation of a Liouville string as a critical equilibrium 
string living in a space time with one extra dimension. In that case the corresponding 
/3-functions of the Liouville-dressed theory would satisfy classical equations of motion. As 
mentioned above, in our approach the conditions describing the restoration of conformal 
invariance by means of Liouville dressing are not to be interpreted as classical equations 
of motion of a string living in a space-time with one extra target dimension. Thus our 
analysis below should be distinguished from previous analysis on Liouville cosmology [ 14j| . 



A generic feature of Liouville dynamics, is that in terms of the world-sheet action, the 
normalization of the Liouville kinetic term can always be arranged (by choosing appropriate 
counterterms) to correspond to a target spacetime of Friedman-Robertson- Walker (FRW) 
type; i.e. the time-like metric component (under the assumption that the Liouville mode 
is time) is: 

Goo = -1. (2.7) 

We remind the reader that the Minkowskian signature is due to supercriticality (3Q 2 = 
c — c* > 0) assumption. This will be understood in what follows. 

3. Non-Criticality Induced in Type-0 String Theory 

We will start our discussion with the ten-dimensional type-0 string theory. The effective 
ten-dimensional target space action of the type-0 String has the form, to 0{a!) in the Regge 



slope a' |jj: 

S = Jd 10 xV^G{e~ 2 *(R + 4(d M $) 2 -^(d M T) 2 -±m 2 T 2 

1 1 T 2 

-y^ H MNp) ~ ^ +T + —)\Fmnpxt\ 2 \ (3.1) 

where capital Greek letters denote ten-dimensional indices, $ is the dilaton, H MNP denotes 
the field strength of the antisymmetric tensor field, which we shall ignore in the present 
work, T is the tachyon field of mass m 2 < 0. In our analysis we have ignored higher-order 
terms in the tachyon potential. The quantity J-'mnpt.t denotes the appropriate five-form 
of type-0 string theory, which couples to the tachyon field in the Ramond-Ramond (RR) 
sector via the function f(T) = 1 + T + \T 2 . 

From ( |3.1| ) one sees easily the important role of the five-form T in stabilizing the ground 
state. Due to its special coupling with the quadratic T 2 term in RR sector of the theory, 
it yields an effective mass term for the tachyon which is positive, despite the originally 
negative m 2 contribution |lj. 

To go off criticality we add the following term to the action (|3.1|) 



J d^V^Ge-^Q 2 (3.2) 



The sign of Q 2 is positive if one assumes supercriticality of the string, which we shall 
do here |I5|, |5], Q. There may be various reasons for having such a term in the effective 
action. One concrete example of induced super-criticality in brane-world scenaria is the 
presence of impulse action (recoil) on the D3-world brane as a result of either quantum 
fluctuations or scattering with other heavy defects or (a macroscopic number of) closed 
string states 0. As argued in those works the impulse deformations of the pertinent 
o"-model are relevant from a world-sheet renormalization-group view point, and as such 
require Liouville dressing ||. The Liouville mode may be identified with the target time 
in a way consistent with the restoration of conformal invariance. In general, there may be 
other microscopic reasons for the induction of non-criticality in string theory, all related 
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to some sort of non-equilibrium processes, which are expected to play a dominant role in 
early Universe Cosmologies. For the purposes of the present article, therefore, we shall not 
specify further the precise microscopic origin of such non-criticality, apart from associating 
it generically with fluctuations of the brane worlds. In this sense we may assume Q to be a 
rather 'phenomenological' parameter, to be determined by consistency with the conformal 
invariance conditions of the Liouville-dressed || a-model. For our purposes we shall restrict 
ourselves to 0(a!) string effective actions. 

In general, Q 2 depends on the a-model backgrounds fields, being the analogue of 
Zamolodchikov's C-function [PJ| . Assuming renormalizability, implies that any explicit 
time-dependence of Q will be absorbed in renormalized couplings g 1 : 

±Q=^Q + P I d I Q = (3.3) 

However, the true set of non-marginal couplings in string theory is actually infinite. More- 
over, in our case, there are non conformal contributions from world-sheet boundaries as 
well, expressing, for instance, recoil of the (quantum) fluctuating D3 brane world. 

Formally, such fluctuations may be approximated semi-classically by impulse-type log- 
arithmic conformal deformations at a a-model level 0, which have recently been shown to 
exist also in theories with world-sheet fermions M , such as the type-0 string theories under 
consideration ffl. Such contributions can only be computed approximately at present, for 
weakly coupled strings 0, §, |16| , and thus are not taken into account explicitly in the effec- 
tive 'phenomenological' approach adopted here, where only the bulk (closed-string sector) 
contributions will be considered in a detailed manner. We denote such contributions by an 
upper case small Latin index {g 1 }, in contrast with the upper case capital index g 1 which 
denotes the complete set of non-conformal deformations. 

In view of ( |3.3| ) this implies that in our case we should consider Q(t,g l ) as depending 
explicitly on time (renormalization-group scale t), in addition to its implicit dependence 
through the renormalized couplings g l . This explicit time dependence simply expresses the 
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contributions to Q from the rest of the non-conformal modes. Its form, then, as a function 
of time t, should be determined by demanding consistency with the generalized conformal- 



invariance conditions (|2.5| ) obtained for the problem at hand from the action ( p.l|) . This is 
the approach we shall assume from now on. 

To this end, we need first to consider the dimensional reduction of the ten-dimensional 
action to the four-dimensional space-time on the brane. This procedure is achieved by 
assuming the following ansatz for the ten-dimensional metric: 



G 



MN 



9 (4) 





e 2ai 
e 2a2 l 



5x5 



(3.4) 



where lower-case Greek indices are four-dimensional space time indices, and I^ x ^ denotes 
the 5x5 unit matrix. We have chosen two different scales for internal space. The field a\ 
sets the scale of the fifth dimension, while CT2 parametrize a flat five dimensional space. In 
the context of homogeneous cosmological models, we are dealing with here, the fields gjfj 
are of Robertson- Walker type, and <Tj, i — 1, 2 are assumed to depend on the time t only. 

By varying the effective action with respect to the five-form one obtains the correspond- 
ing /^-function: 



ft 



M 



[l+T 



T 2 



)T 



MNPET 



(3.5) 



It is important to notice that in our approach we consider as a source of non-criticality for 
the underlying string theory only quantum fluctuations of the D3 brane resulting in the 
non- vanishing of the renormalization group /3-functions of the string multiplet, i.e. graviton, 
dilaton and tachyon fields. The deformation corresponding to the five-form |l] is assumed 
conformal. This is associated with that fact that this field is related to the RR electric flux 
of the D-branes, and as such is not affected by the (recoil) quantum fluctuations. Moreover, 
for our purposes in this work, which is to demonstrate an inflationary phase and graceful 
exit from it, the essential role is played by the string multiplet fields, with the five-form 
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flux field serving the sole purpose of tachyon stabilization. This implies the vanishing of 



the left-hand side of (|3.5|): 

T 2 



= V 



M 



;i + r + — )/• 



MNPT,T 



(3.6) 



Moreover, for simplicity we shall assume a configuration for which the five-form acquires 
non- vanishing values only along one extra dimension: 

FMNPUT = £0af3-y8 /5W (3-7) 

where the lower-case Greek letters denote four-dimensional indices. 

The ten-dimensional Weyl-anomaly coefficients (generalized /3-functions |)17[| ) for the 
cr-model corresponding to propagation in a type-0 string background read, to 0(a'): 



1 



v2 in-r tKAPZ 



7} e f(T) Gmn{Fk\PY,t) — lO^Mi^APE-^Af 

~f = {-V 2 T + 2T ;M $ M + m 2 T} + 2e 2 *f(r)(^MVPET) 2 (3.8) 

Due to the renormalizability of the (non-critical) a-model, there is an additional equa- 



tion |L8| which should supplement equation (|2.5| ), the Curci-Paffutti relation, which relates 
the dilaton /3-function, and hence the effective running central charge of the theory, with 
the rest of the /3-functions: 

- V M /3* + 2G NP e 2 *V N (e- 2 *pf IP ) + T, M ~f = (3.9) 

Although this equation holds formally in the flat world sheet case, however in our frame- 
work it should also hold for the (3 % {X) functions, i.e. the flat world sheet /5-functions upon 
the substitution of the a-model couplings with the Liouville dressed ones. It provides a 
highly non-trivial constraint Q, which should be respected by the process of identifying the 
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Liouville (world-sheet) renormalization scale with the target time ||. As emphasized by 
Tseytlin |L7| this equation expresses, under the off-shell equivalence of the Weyl-anomaly 
coefficients f3 % with the field variations obtained from the target space string effective ac- 
tion, the invariance of the latter under general coordinate transformations. As such, the 
extension of this consistency condition to our Liouville string case, with the Liouville mode 
playing the role of the target time, is evident. Notice that, due to the criticality condition 
( J3l)l) , the five- form /3-function does not appear in ( [31)1) . 

An important comment is now in order concerning the role of the dilaton field in non- 
critical string theories. In general, it is not clear that the dilaton can be treated as one of 
the other couplings g % of the deformed a-model. In the Appendix we specify the precise 
conditions under which this can indeed be made possible, thereby allowing the application 
of fl2.5p to the dilaton coupling as well. 

Upon considering the fields to be time dependent only, i.e. considering spherically- 
symmetric homogeneous backgrounds, restricting ourselves to the compactification (|3.4j) , 
and assuming a Robert son- Walker form of the four-dimensional metric, with scale fac- 
tor a(t), the generalized conformal invariance conditions (|2.5|) , using ( |3.8| ),( |3~9D imply the 
following equations: 

-3- +*i + 5a 2 - 2$ + &{ + 5d 2 + -T 2 + e- 2CT1+2 */ 2 /(T) = , 
da 



+ aa (2Q +& x + 5a 2 - 2$) + e- 2 ^ 2 * f 2 f{T)a 2 = , 
3a! + 5d 2 + 3-cti + 2Qcxi + 5dia 2 - 2ai$ + e" 2,7l+2 */ 2 f(T) = 



3a 2 + 9d 2 + 3-a- 2 + 2gd 2 + a-ia-2 - 2d 2 <i> - e - 2CT1+2 */ 2 f(T) = 

CI 

if + 3-T + Q f + &xf + 5a 2 T - 2T$ + 
a 



m 2 T - 4e- 2CT1+2 */ 2 /'(T) = , 



a „a 2 



$ + g$ + 6- + 6— + 



a a 2 
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— <7i 0"i — 5(72 — 15 — 0"2 — 0" 1 — 15<7 2 — 5(71(72 — 

a a 



2 <j> 2 + 2$ + 6-$ + 2cr!$ + 10d 2 $ 
a 

C 5 = e-°^f{T)h , 



^T 2 + ^m 2 T 2 + Q 2 = 
44 



$(3) + g$ + q$ + 12 _^ (oo + d 2 + gad)- T(T + QT) + 

4di(ai + 2d 2 + Qo-i) + 20a2(<72 + 2d 2 + Q& 2 ) = (3.10) 

where f'(T) denotes functional differentiation with respect to the field T, the overdot 
denotes time derivative, and $ < - 3 * ) denotes triple time derivative. 

Note that the above equations have been derived from a ten-dimensional action. An 
equivalent set of equations (in fact at most linear combinations) come out from the cor- 
responding four- dimensional action after dimensional reduction. Of course this reduction 
leads to the string frame. We may turn to the Einstein frame through the transformation 

g E = e -2*+-i+5^ (3 U ) 

In this frame the line element is 

ds 2 E = -e- 2 * + ^ +5 ^dt 2 + a 2 (t)e- 2 *+ CT1+5CT2 (dr 2 + r 2 dQ 2 ) (3.12) 

Therefore to discuss cosmological evolution we have to pass to the cosmological time defined 

by 

Then the line element becomes 

ds% = -dt\ + a 2 (t E )(dr 2 + r 2 dVt 2 ) (3.14) 

with a'(t E ) = e-* +zl± ^a(t(t E )). 

We should also remark that we have adopted the usual Kaluza-Klein reduction without 
considering the extra dimensions compactified. Nevertheless the equations we are interested 
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in and the results we will discuss in the following section are not modified if we had 
considered a compact six- dimensional space instead. In that case e 2cri and e 2<J2 should 
correspond to radii of the compact space. Note also that in the string frame we have the 
exponential e - 2 *+°"i+ 5o "2 j ns tead of e -2 *, since we allow time dependence of the volume of 
the compact space. 

4. Solution of the System of Differential Equations 

The problem of initial conditions of the Universe is a very complicated problem in 
Cosmology, and the present four- dimensional model is no exception. To bypass this problem 
we adopt here a different strategy. We demand that the solution approach asymptotically 
the Minkowski space time (in the sigma- model framework), with a linear dilaton. This 



space time is well-studied from the point of view of non-critical string theory JT5[, and 
provides an example of a linearly expanding, non accelerating Universe in Einstein frame Q 
From the point of view of a local field theory, the fact that the acceleration of the type- 
0-string Universe is not eternal is essential for a consistent quantization, and definition of 
asymptotic states, and hence a scattering matrix |2U], |2T], [22] . 

We therefore assume that the universe of a type-0 non-critical string theory is described 
by a configuration, which, in the u-model frame, asymptotically (i.e. for large times) 
approaches a Minkowski spacetime with linear dilaton and constant values of U\ and o"2, 
and examine whether such an assumption leads to a consistent (and continuous) solution 
of the system of equations ( |3.10|) for all times. For this purpose we separate the fields in 



their asymptotic values plus fields which tend asymptotically to zero. Substituting these 



1 Concerning the model of [|5| there has been some discussion in the literature [|1| as to whether this 
type of universe is really (physically) expanding or is actually static, in case one accepts the point of view 
that the length measurements are made by string rods. In our context there is no such ambiguity, because, 
as we shall see below, this type of Universe is the final stage of the Liouville-time evolution, and for us the 
Einstein frame, in which one recovers the linearly expanding Universe of JL5|, is the physical frame. 
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fields back to the equations we let the system evolve in time backwards. 

We now describe the various steps of this procedure in detail. The separation of the 
various fields in their asymptotic values plus fields which tend asymptotically to zero, is 
performed as follows: 

$(£) = f + fit + h(t) 
a{t) = a Q e m 
cri{t) = S l + Si(t) 

02(*) = S 02 + S 2 (t) 

Q{t) = q + qi(t) 
T{t) = c + r (t) 
F 5 (t) = e*f 6 (t) (4.1) 



where the constants go and /i are related through the relation q = — & \l + Vl7}. This 
relation results from the requirement that the dilaton equation is satisfied. The fields 
I h(t), b(t), si(t), S2(t), qi(t), T (t), f 5 (t) \ vanish asymptotically. We want gravity to be weak 
asymptotically so we choose /i to be negative. Due to the bare tachyonic mass, m° < 0, 
Co must be chosen to vanish for the tachyon equation to be satisfied asymptotically. Note 
that the above asymptotic conditions in the Einstein frame cosmology correspond to a 
linear expanding universe |Tj| . With this choice the system of the equations ( |3.10| ) become 
quasi-linear and, in order to reduce the order of the system, we define 

h(t) = h (t) 

h{t) = hx{t) 

b{t) = b {t) 

b{t) = h(t) 
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si{t) = s w (t) 

s'i(t) = sn(t) 

s 2 (t) = s 20 (t) 

s 2 (t) = s 21 (t) 

T(t) = %{t) 

t{t) = T x (t) (4.2) 

So the system of equations ( ft.lOQ take the form 

S=Ax + F{x) (4.3) 

where x = {h (t), h^t), q^t), b (t), h(t), s 10 (t), s n (t), s 20 (t), s 2 i(t), %{t), T x (t), / 5 (t)| , with 
_L denoting transpose, A is the constant matrix determining the linear part of the system 
and F(x) gives the nonlinear terms. 
The differential system (|4.3|) can be written formally as an integral system 

x{t) = x(t ) + Y(t) /* dsY-\s)F [x(s)] , (4.4) 

J to 

where the matrix Y is a fundamental solution of the linear system satisfying the equation 

Y = AY (4.5) 

The full solution of this integral system can be given in an iterative form: 

%+!)(*) = %)(*) + Y(t) f dsY-^F \x {n) (s)} (4.6) 



The starting point of the iteration procedure is the solution of the linear system with 
the correct asymptotic behaviour. If we insert the fields ( |4.1| ) and ( |4.2| ) into the system of 
equations (|3.10|) and keep only the linear part we get the following general solution of the 



system 



h{t) = Cho + Chtf + \ (C Jl + 5C , S2l) e*<*-«* 

4 (/i - go) 



where 
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2\pfo-soi-5so2+2fit 



Qi(t) 



b(t) 

s 2 (t) 



3 Ch 2(h _ qo)t C 5 Cf 5 (Qtf + 9f iqo + q 2 )e 

4 (A -go) 8/ igo (2/i + go) 

3Cb 1 (2f 1 -q ) 2(h _ qn)t 
2/i 

(Cs n + 5(7320(2/! + go) | (/l - go)i 
6/1 

C 5 Cf 5 (6ff + 9/ igo + g 2 ) /n _ Sm _ 5 , n2+2/lt 
8A 2 go 

(^A 1 ^1 p 2(/i-g )i _ C5C/5 / -aoi-5a 3+2/il 

G00+ 2(/ 1 -g ) 8/ igo 

C Sl0 + _i B ££iL_ e |(/i-9o)t _ ^C/s _ c / - B01 -5802+2/^ 



2(/i - go) 



8/1(2/1 + go) 



C 20 + 3gg 21 c f(/ 1 - qo )t + C5C/5 e f -a 01 -5ao2+2fit 



2(/i - go) 



8/1(2/1 + g ) 



ht 



Che 



Cl e M + C 2 e Bt + 



C 5 Cf, 



8/2 + 2m 2 + 4/igo 



e /o-soi-5so2+2/i< 



(4.7) 



.4 
B 



(2/i-g )V(2/i-go) 2 -8m 5 



(2/i + go)V(2/i-go) 2 -8m 2 



<0 

>o, 



(U 



and jClio, CTii, Cgi, Cb , Cb\, Csio, Csn, CS20, CS22, Ci, C2, C/s| are constants of integra- 
tion. To have the right asymptotic behaviour we choose lCho,Chi,Cbo,Cqx,Csio,Cs2o\ 
to be zero. The constant C 2 corresponds to a positive eigenvalue of the matrix A, so we set 
it also zero. Since the constant matrix A has negative eigenvalues, well known theorems 
from qualitative theory of quasi-linear systems, guarantee the asymptotic stability of the 
iterative solution 



Now at the first step in the iteration the solution is given by: 



x(t) = X( 0) (t) +Y(t) / dsY- 1 (s)F Y(s)C 



to 



(4.9) 



where X(p)(t) is given by, 

x (0) (t)=Y(t)C (4.10) 

and C is a constant vector C = \Ch , Chi, Cqi, Cb , Cbi, Cs w , Csu, Cs 2 q, Cs 2 2, C\, C 2 , C/ 5 \ 
with Cho = Chi = Cq± = Cbo = Csiq = CS20 = Ci = for reasons we explained above. 

The first step of the iteration can be integrated analytically and has the right asymptotic 
behaviour as expected. Nevertheless due to the complexity of the expressions and the fact 
that they do not give us the full solution of the problem, we proceed as follows. We take 
from this step the value of the fields at late times and we use them as initial conditions for 
the numerical solution of the system. 
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Figure 1: The evolution of the scale factor in time. 

We can now discuss our results. As we can see from Fig.l-Fig.3, the cosmological evolu- 
tion of our universe passes through the following phases. At very early times, the universe 
starts from the initial singularity. Then enters a phase where the physical dimensions are 
formed. At this stage <Ji, a 2 and a are comparable in magnitude. The phase of inflation 
follows, during which the scale factor grows exponentially (for a short time though), while 
the internal space contracts with very negative values of a\ and <j 2 - Finally, the universe 
enters a phase, where its expansion slows down until it reaches the asymptotically flat space 
(in the a- model frame). The internal space continues to contract but with very slow rate 
(with &i and &2 positive) until it reaches a constant value. As we can see from Fig. 2 and 
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Fig.3, o\ and 02 scale differently. The field <72 very soon freezes to a value much lower than 
the value of a\, indicating the the fifth dimension can be much larger than the other five 
dimensions. 

si 
-0.6' 




Figure 2: The behaviour of the field representing the evolution of the fifth dimension. 
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Figure 3: The behaviour of the field o<i expressing the evolution of the five internal dimen- 
sions. The difference of the scales in the decomposition of the extra dimensions is clearly 



seen. 



For the dilaton field (Fig. 4), we observe that at the singularity the dilaton field is 
infinite, indicating that the gravity is very strong. Then at the second phase of evolution, 
the strength of gravity is weakened because the dilaton field drops linearly. This weakening 
of the gravitational interactions continues during inflation, and finally at the exit, the 
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dilaton field continues to drop linearly. The Tachyon field (Fig. 5) falls continuously during 
the evolution until it becomes zero. The RR-field (Fig. 6), on the other hand, grows 
from zero value at the initial singularity until it reaches a constant value at the exit of 
inflation. Moreover, the behaviour of the central-charge deficit Q 2 is depicted in figure 7. 
The latter is initially decreasing, until it becomes zero, and then oscillates before reaching 
an asymptotic constant positive value Q 2 — > const ^ 0, as t —* oo. A similar behaviour is 
observed in the two dimensional case ||. Notice that the oscillations in the central charge 
deficit are compatible with the fact that the Liouville mode is not a unitary a-model field, 
as a result of its negative signature. Hence, Zamolodchikov's c-theorem does not strictly 
apply, although overall there is a reduction in the central charge, since the theory starts 
with a very large value of Q 2 (formally infinite near the initial singularity), which eventually 
asymptotes to a finite positive value. 




Figure 4: The behaviour of the dilaton field. 



We can analytically support the numerical solution we found looking more carefully at 
the equations satisfied by the fields. Assuming that asymptotically the fields o\ and a-i are 
constant, using the second and seventh equation from the system of equations ( |3.10| ), the 
field b takes the form of a Riccati equation, which is the usual case in cosmological models 
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Figure 5: The behaviour of the tachyon field. 




Figure 6: The evolution of the field f& denoting the behaviour of the R-R flux. 



based on type-0 string theory H, 



C, 



b + b(2Q - 2$ + Ji + 5<r 2 ) +b 2 = —±e 



-<Tl+5(T2+* 



(4.11) 



For times nearly after the initial singularity and during inflation, equation ( 4.11|) is approx- 
imated by 

b + 2b(Q(t) - $) + b 2 = (4.12) 

This equation has the solution 



b = c\ + t + ln(l — e 



-cat\ 



(4.13) 
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Figure 7: The behaviour of the central charge deficit Q 2 . This quantity, after some oscil- 
latory behaviour near its zero value, asymptotes for large times to a constant value (not 
depicted clearly in the figure due to numerical limitations). 

where C\,C2 are constants. This solution shows that the scale factor has an inflationary 
phase and for later times the logarithmic term takes over and slows down the expansion of 
the universe. 

For large times when the internal dimensions have been stabilized in constant values, 
the equation ( |4.11| ) becomes 



b + 2b(q - A) + b 2 = ce 2ht (4.14) 



The exact solution of this equation respecting the required asymptotic behaviour is 



b(t) 



In [0F1 (|, ^e 2 ^) + c 3 e 2 ^ V x ( 2 - |, ^A] (4-15) 



where C3 is a constant, and o-^i is the Confluent Hypergeo metric Function. If we expand 
this solution in series, we get the asymptotic behaviour of b(t) which characterizes the 
inflationary exit period. The quantity Q(t) — $ is crucial for the existence of an inflationary 
phase and then exit from it. We note that in two target dimensions, it is the same quantity 
that controls the exit from inflation [HI. 
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We now make some brief remarks on the physical significance of the asymptotic be- 
haviour of our solution (for large times t —>■ oo). A detailed discussion will appear in a 
forthcoming publication |24j]. By inspecting figure 1, or analytically from Q4.15 ), we ob- 



serve that the acceleration of the universe vanishes asymptotically in time (t — > oo) and 
one recovers a non-accelerating Minkowski universe with a linear dilaton in the cr-model 
(string) frame |15|]: da = b + (b) 2 — > 0, for t — > oo, with a(t) — > a . 

Such configurations imply the absence of a cosmological horizons, given that the relevant 
quantity diverges, 5 oc J^ dt/a(t) — > oo |. This is a welcome fact from the point of view of 
low-energy field theory, in the sense that in such not-eternally accelerating Universes one is 
able to define asymptotic states and thus a proper S matrix |pCfl . Out of interest we notice 
that asymptotic situations of this kind appear to be generic in cosmologies with fluctuating 
D-brane defects |7|]. 

Another issue that arises in this scenario is the asymptotic value of the vacuum energy. 



In the Einstein frame ( |3.13|) , the latter is given by 



where we took into account that asymptotically in time the fields a±, a 2 become constants, 
the dilaton configuration is linear and negative (see figure 4), while the central-charge 
deficit Q 2 -^ const > 0. In our solution, therefore, we observe that the four- dimensional 
vacuum energy is positive (de Sitter type) and tends to zero asymptotically. When we 
express the result (|4.16| ) in terms of the Robert son- Walker time t^ ( |3.14|) we observe that 



the four- dimensional vacuum energy relaxes to zero as: 



A = (const)— (4.17) 



ti 



2 Notice that a divergent result for S is, of course, also obtained in the Einstein frame ( 3.13|) , ( |3.14| ) 



when we express the integrand in terms of the Robertson- Walker time £#. In that case the scale factor is 
expanding linearly a'(ts) oc i_e, while the dilaton varies logarithmically with the time £#. 
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which is compatible with recent experimental evidence, and is obtained in many quintessence 
models. In our approach the quintessence field is, in a sense, the dilaton |13|, ^TJ. We shall 
present a more detailed account of these features in a forthcoming publication [Q . 

We next remark that in the context of a critical string theory, with Q = 0, the dilaton 
equation in ( |3.10| ) implies that the dilaton field $ becomes asymptotically constant and, 



hence, equation ( [4. 14] ) becomes 

b + b 2 = c (4.18) 

This has as solution 

b(t) = ci + In coshU/cb - y/cc 2 \ (4.19) 

where c\ and c 2 are constants. Then the scale factor has an exponential growth but it will 
expand for ever. Such a case is inconsistent from the point of view of defining asymptotic 
states and an S'-matrix. Because the case Q = 0, is supposed to correspond to critical string 
theory, which must be a theory of an S'-matrix, we may conclude that such a situation will 
not arise as a consistent string theory solution. 

All the analysis we have obtained so far pertains to the string (cr-model) frame. As 
already mentioned, one gets similar results in the Einstein frame, using ( fj . 1 3| ) and ( |3.14| ). 

Another consistency check of the approach is provided by the F-field of the RR-flux. 
Because of equation ( |3.6|) , the field F 5 (t) has the exact solution 

F 5 (t) = e*h = ~r^ 7^e flt+h{t) (4-20) 

Stability of the internal space and consistency of our scheme requires the constant fo to have 
a large negative value. This means that gravity is asymptotically very weak, as expected. 

5. Discussion 

In the present article we have presented a cosmological model based on a type-0 string 
theory. The type-0 string theory is rich in its content. In addition to the graviton and 
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dilaton fields, it includes also a tachyon field which couples to an RR five-form field. To 
avoid tachyon instabilities the tachyon field has to take on values at the minimum of its 
potential. Considering the ten-dimensional action of the type-0 field theory one can get 
at the conformal point the /3-functions of the theory, and from them, by reduction to 
four-dimensions, and assuming that all fields are time dependent, one obtains an effective 
four-dimensional theory. It is found that this theory, in the context of a Robertson- Walker 
background, has an inflationary phase but there is no smooth exit from this phase. 

We believe that the central issue of inflationary cosmology is not how one can get an 



inflationary phase, but how one can exit from it. In view of the recent discussion [20| of 



the exit from inflation in string theory, and in general on the consistent quantization of 



de-Sitter Universes |22| , we think that it is useful to propose new mechanisms to exit from 
inflation. In our previous work |JJ, we had proposed a mechanism based on non-critical 
strings, to exit from inflation. We had considered a two dimensional model with a tachyon 
field. The presence of the tachyon field was crucial for an inflationary phase to go smoothly 
to a Robertson- Walker phase. Notably, in two-dimensions the tachyon field is an excitation 
of the 'massless' matter multiplet, and is not associated with an instability of the vacuum. 

In the present work we have applied this mechanism to a realistic four-dimensional 
model of type string theory, where again the tachyon is stabilized by a RR flux field. We 
modified the /3-functions of the ten-dimensional theory, assuming that the string theory is 
non-critical, which in such theories may be caused by fluctuations of the 3-brane worlds. 
This hypothesis introduces new terms in the /3-functions of the theory. Physically these 
terms express the fact that our non-critical string theory is performing small oscillations 
around the conformal point. 

The modified /5-functions supplemented with the Curci-Paffuti equation implied by 
world-sheet renormalizability, have been reduced to four- dimensions. Assuming a homo- 
geneous and spherically symmetric background, we have solved numerically the resulting 
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equations. Our solution demonstrates that the scale factor of the Universe, after the ini- 
tial singularity, enters a short inflationary phase and then, in a smooth way, goes onto a 
flat Minkowski spacetime with a linear dilaton for t — > oo. The fields <j\ and a 2 which 
parametrise the internal space have an interesting behaviour. The field <Ti, which sets the 
scale of the fifth dimension, during inflation contracts until it reaches a constant value. 
After inflation, it maintains this value, until the universe evolves to an asymptotically 
flat spacetime (in the string frame, or a linearly expanding, non-accelerating Universe in 
the Einstein frame). The field a 2 which parametrises the conformally flat five-dimensional 
space freezes to a constant value which is much smaller than that of the fifth dimension. 
Thus we see that, in our model, a cosmological evolution may lead to different scales for 
the extra dimensions. 

It is important to notice that the contraction of the extra dimensions is due to the 
fact that gravity is very weak asymptotically in our model. It is also worthy of stressing 
that the possibility of having one large extra dimension is achieved upon the choice of the 
RR flux along that direction. Hence, such a possibility arises in our model from string 
non-perturbative effects. A phenomenologically important feature of the model is that the 
vacuum energy, determined by the central-charge deficit, relaxes to zero asymptotically in 
a way which is reminiscent of quintessence models, with the role of the quintessence field 
played by the dilaton. 

The type-0 string Universe does not accelerate eternally, thus avoiding the problem 
associated with the presence of cosmological horizons, as far as the definition of proper 
asymptotic states, and hence an S-matrix, is concerned. A detailed study of such issues, as 
well as of re-heating and other physical features of inflationary models, will be the subject 
of forthcoming work. 

We close our discussion with an important comment concerning the initial singularity 
of our solution. The singularity is a general feature of the equations of the form (|3.10p . 
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Removing the singularity is probably a matter of a full quantum description of the theory, 
which at present is not available. We note at this stage, however, the possibility of deriving 
smooth cosmological solutions in string theory, without initial singularities, by including 
in the action higher curvature terms (e.g. quadratic of Gauss-Bonnet type f2~5fl), which are 
part of the quantum corrections, in the sense of being generated by including string-loop 
corrections in the effective action. We plan to return to a systematic study of the effects 
of such higher-order terms, within our framework, in due course. 
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Appendix 

Dilaton Equations in Liouville Strings 

In this Appendix we explain in some detail why the dilaton field obeys an equation 
of the form ( |2.5|) , in a similar manner with the rest of the background fields/couplings of 
the deformed cr-model. At a a-model level, the dilaton field couples to the world-sheet 
curvature in the same way as the Liouville mode, which implies that in general the dilaton 
of the Liouville-dressed theory will receive contributions from the Liouville mode p. The 
Liouville mode p is viewed as an extra target space dimension, which has Minkowskian 



signature in supercritical theories fl5|| . In our approach, we first treat this extra time 
dimension as a second time, which, however, is eventually going to be identified with the 
existing time-coordinate t of the non-critical string. 

In this sense, the restoration of conformal invariance implied by the Liouville dressing H 
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means that in D + 1 dimensions, with D the original dimensionality of the non-critical 
theory, one would have a conformal string. In terms of the /?* function we would have: 
P$(D+i) _ q(^ — 10 = 0, where C tot is the total central charge of the Liouville-dressed 
theory, and the superscript D + 1 denotes (D + l)-dimensional quantities. To order 0(a') 
in a Regge slope a' expansion, one has for a <i-dimensional target space: 



Ctot = d- -a 



R - 4(V$) 2 + 4V 2 $ + j{S/Tf + \m 2 T 2 



(5.1; 



In general, for non-critical strings with central-charge deficit Q 2 , one should replace d by 
3(Q 2 + d*), where d* is the critical value of the central charge for the conformal (fixed-point) 
theory (d* = 10 for superstrings we are dealing with here). For convenience we set from 
now on a' = 2, which is the normalization used in (|3.8f) and throughout this work. 

One now recalls the definition /3*( D+1 ) = / 3*( D + 1 ) — G pu f3^ D+1 \ where p,u span the 
D + 1-dimensional spacetime, including the Liouville dimension p. In Liouville theory, the 
metric G pp = — 1, G p m = 0, and hence j3^ = = fl? M , which implies that the Liouville 
renormalization does not introduce any extra contributions to the trace of the graviton (3 
function, and therefore to the space-time curvature terms in the expression for /?* (|3.8|) . 
This leaves one with the following condition, expressing restoration of conformal invariance: 



AG PP 



(9 P $)(9 P $) - d p d p ® - l -G MN {d p G MN ){d p §) + l(d p T)(d p T) 



+ pHD) = Q ( 52 ) 



where (3®^ refers to the D-dimensional parts of the dilaton Weyl anomaly coefficient 
which notably include the dilaton field. 

As we shall discuss below, the terms \G MN d p GMNd p <& may be neglected. There are 
two ways in which this can be justified. First, in our approach we treat the Liouville 
mode p as a world-sheet renormalization-group scale. Hence terms like 8 p Gmn = Pmn 
and d p Q = (3 coincide with the world-sheet graviton and dilaton /5-functions. Thus, such 
terms are of higher order, as compared with the rest of the terms in (|5.2j ), if one works in 
the neighbourhood of a fixed point, which we assume here. Alternatively, such quadratic in 
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the /3-functions terms can be eliminated by an appropriate renormalization-group scheme 
choice fl26f , as we shall discuss later. 



We now assume that the Liouville mode p contributions to $(p, t, . . .) can be represented 
by 

${p,t i ...) = -<yQ(p)p + <p{p i t,...) (5.3) 

where the . . . denote possible dependence on spatial target coordinates, 7 is a numerical 
constant, to be determined below, and the term proportional to Qp is motivated by the 
Liouville dynamics || |5[. Under this splitting, and taking into account that Q does not 
depend on any other coordinate except the local renormalization-group scale p || , one may 
write (|5.2j ) for the supercritical string case G pp — — 1: 

ACp - 4 7 2 Q 2 + 8 7 Q0 = -0* + 0[(d p <5>) 2 , (d p T) 2 , d p Q, (d p G)(d p $)] (5.4) 



where d p Q oc (f3 1 ) 2 by virtue of the Zamolodchikov c-theorem fLl] , extended to the local 
scale p ||. 

Then, terms of the form 0[(d p (p) 2 , (<9 p T) 2 , d p Q, (d p G)(d p &)] can be removed by field 
redefinitions (appropriate renormalization-group scheme choice). This is possible due to 
the gradient flow property (2J5) of the j3- functions [Pql , which can be shown [§ to charac- 



terize the Liouville problem upon considering || the Liouville mode as a local scale on the 
world-sheet, appropriate for curved- (two-dimensional)-space renormalization |12|]. Indeed, 



consider an infinitesimal field-redefinition 5g l (or equivalently 5X l for the Liouville-dressed 
couplings), where now the set g l includes the dilaton field. This expresses a particular 
choice of (a class) of renormalization group schemes. Under this redefinition, the central 
charge deficit Q 2 , which is assumed to satisfy the gradient flow property di{Q 2 ) = Gijft , 



changes as 



5(Q 2 ) = Sg%(Q 2 ) = SgVijP 3 (5.5) 



It is then clear that by choosing appropriately Sg l oc /3 l , one may absorb terms of the 
form 0[(d p ip) 2 , (d p T) 2 ,d p Q, (d p G)(d p $)] on the right-hand-side of Q f26|. Moreover, as 
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already mentioned, such terms may also be neglected due the fact that, in the models of 
'recoil' fluctuations of D3-brane worlds we are dealing with here |7|, ||, the induced non- 
criticality is marginal, thereby forcing us to work in the neighbourhood of a world-sheet 
renormalization-group fixed point. Hence, terms quadratic in the j3- functions, and thus in 
d p X l , are subdominant compared to terms linear in such quantities. 

With the above in mind, we now observe that, upon identifying p = 2t, where t is the 
target time, fixing 7 = 1/4, and re-expressing the terms involving tp in terms of $, we may 
rewrite (|5.4|) in the form (|2.5|) , up to terms (denoted by . . .), which may either be removed 
by a renormalization-scheme choice or are negligible in magnitude for an expansion in the 
neighbourhood of a fixed point we are assuming throughout: 

= (p + Qp + 0* - 1q 2 = $ + Q$ + P* + . . . = (5.6) 

The overdot denotes differentiation with respect to time t. Notice that the identification of 
the local world-sheet renormalization-group scale p with the target time may be understood 
as a renormalization-group scheme choice ||. 

Thus, in this way, the dilaton equation (|5.6|) is now expressed as an equation of the form 
( 2-5 ). Therefore, by an appropriate scheme choice, one was able use the equations ( |2.5|) for 



all cr-model couplings, including the dilaton, in a way consistent with both the restoration 
of conformal invariance by the Liouville mode, and the identification of the latter with the 
target time. 
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